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SUMMARY 


More stringent noise regulations for commercial jet aircraft have prompted continuing research into 
methods for the design of improved engine noise suppression systems. This report presents the 
results of a research program that was undertaken to analytically and experimentally evaluate the 
use of phased multisegment linings to provide more efficient attenuation than that provided by 
current generation linens. Analytical methods were developed to take advantage of relative 
placement of lining segments so as to reflect and spatially redistribute the acoustic energy incident 
on the lining. These multisegment linings were shown to outperform single-segment optimum liners. 
The results of this investigation are immediately applicable to the design of improved aircraft 
exhaust duct noise suppressors. Further investigations appear to be necessary before phased liner 
techniques can definitely be shown to be adaptable to inlet liner design. The overall Investigation 
involved analytical studies conducted to identify the principal mechanisms by which phased linings 
achieve improved attenuations, and the development of methods to utilize these mechanisms in the 
design of improved liners which were built and then tested in a laboratory facility. The test data 
were subsequently analyzed both to substantiate and to augment the analytical techniques. 

The analyses were based on a mathematical model for rectangular ducts with uniform mean flow. 
Segmented acoustic fields were represented by duct eigenfunction expansions and mode-matching 
was used to ensure continuity of the total field. Parametric studies were performed to identify 
attenuation mechanisms and define preliminary liner configurations. An optimization procedure 
was used to determine optimum liner impedance values for a given total lining length, Mach 
number, and incident modal distribution. Optimal segmented liners were developed and it was 
shown that, provided the sound source is well-defined and flow environment is known, conventional 
infinite duct optimum attenuation rates can be exceeded. The primary suppression mechanism of 
the multisegment liners was found to be "modal conditioning." {The incident sound is redistributed 
by the front liners into modes that are more readily absorbed by the remaining lining segments.) 

An experimental program was conducted in a laboratory test facility to confirm the analysis results. 
The measured data are presented in the form of analytical-experimental correlations. Six lining 
configurations of one, two and three segments were designed and tested for two flow conditions, 
M = 0.0 and M = 0.4. The two- and three-segment liners designed for M = 0 produced 4.6 and 6.4 dB 
more PWL attenuation respectively than the single-segment liner for the design frequency. For the 
M = 0.4 tests, the two- and three-segment liners produced 6.0 and 10.0 dB more attenuation than 
the single-segment liner. Results of this research demonstrate that acoustic linings with improved 
suppression effectiveness can be designed for application to future quiet propulsion systems which 
meet more stringent noise goals. 


INTRODUCTION 


Background 

Current and proposed notse regulations for commercial aircraft have prompted continuing research 
and development of improved suppression systems for turbofan engines. Although some reduction 
in turbofan noise has been achieved by suitable engine cycle selection and component design 
acoustic treatment is still required to suppress the internally generated noise. The imposition of 
increasingly stringent noise reduction goals, which must be met with minimum impact on airplane 
performance, requires that acoustic lining design techniques must be continually refined and 
upated. 

Standard uniform infinite duct theories (ref, 1} are based on the assumption that the duct acoustic 
field includes only modes moving away from the source. The acoustic energy in the duct may then 
be uniquely associated with the cut-on or propagating modes. Cremer's theory (ref. 2), which 
involves choosing a liner for which the attenuation rate for the least attenuated mode is optimized, 
has been reinvestigated in detail by Tester (ref, 3) and Zorumski and Mason (ref. 4), 

Rice (ref. 5) has shown that, for an initially plane pressure distribution in an infinite uniform lined 
duct, the optimum liner would produce much less attenuation than that predicted by the Cramer 
theory. This is due to the presence of other modes which are required to produce the initial plane 
pressure wave. The sound power level attenuation rate for the complicated wave system is quite 
different from that of any single mode. The total sound power includes not only direct 
contributions from each mode, but also cross terms contributed by modal interaction. For very long 
ducts the plane wave attenuation per unit length will approach the Cremer value; however, the total 
attenuation over a given length will be different. 

Obvious!/, this method of approach (See also Snow, ref. 6, who investigated other pressure 
distributions in addition to the plane wave) should yield more meaningful results than the Cremer 
theory for the lining lengths typically used in practical applications. However, regardless of which of 
the infinite theories is used, there are still factors related to finite liner length and the impedance 
discontinuity at the boundary between liners which have not been taken into account. 
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Method of Approach 


To evaluate finite liner length effects on attenuation performance, acoustic wave reflection from 
duct nonuniformities should be considered. The computational method generally used for this 
purpose is the mode matching technique, in which the acoustic field is represented by eigenfunction 
expansions in terms of incident and reflected duct modes. The undetermined modal coefficients for 
the field description in any uniform section are determined by requiring continuity of the fields at 
segment junctions. This method has only recently begun to be applied to acoustic investigations 
(refs. 7-12) although it has been often used in electromagnetic applications {refs. 13-15). 

Lansing and Zorumski (ref. 7) were the first to use this theory to investigate the advantages 
afforded by the inclusion of impedance discontinuity effects in the evaluation of liner performance. 
Their preliminary study provided limited insight into the mechanisms responsible for the behavior 
of the linings. Among the limitations of their study were the facts that liner admittance was 
assumed constant with frequency, and purely reactive liners were allowed; realistic liners have 
characteristics which vary with frequency and have finite resistance values. 

The present investigation was conceived as an extension of the work of Lansing and Zorumski (ref. 
7), aimed at accounting for some of the limitations of their investigation. Identification of the 
segmented duct mechanisms and experimental verification of predicted performance were primary 
goals. Realistic liner impedance models were used and very practical constraints on liner design were 
Imposed so that the test panels could be constructed using conventional materials. 

To identify mechanisms and develop liner designs, detailed analytical studies were made. These 
studies involved both extensive parameter variations and numerical optimization runs. Liners chosen 
for testing were built and installed in a grazing flow duct facility. Performance comparisons and 
analytical-experimental correlations were based on comparisons of measured and analytical axial 
distributions of pressure magnitude and phase. The philosophy of the bench-test program was to 
compare optimum single-segment performance with that of the best multisegment liner which could 
be built subject to the constraints of material availability and uncertainty regarding details of sound 
source modal structure. 

The multiparameter optimization technique used in the analytical studies to be reported here was 
the conjugate gradient method. Similar steepest descent algorithms have been used in liner 
investigations by Martenson and Liu (ref. 16) and Wilkinson (ref. 17) and in expansion chamber 
studies by Kessler and Puri (ref. 18). However, the present investigation and concurrent work by H. 
C. Lester and J. W, Posey (private communication) of NASA Langley Research Center appear to be 
the first in which highly accurate mathematical models and large numbers of parameters (lengths 
and impedances of multiple segments) have been considered. 
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Complex modal coefficient for nth right moving mode in Jth segment 
Liner specific admittance 

Value of definite integral defined in Appendix E and G 
Complex modal cv. ‘^icient for nth left moving mode in jth segment 
Value of definite integral defined in Appendix E 
Speed of sound 

Value of definite integral defined in Appendix G 
Denotes differential 

Eigenvalue equation defined in Appendix B 
Denotes base of natural logrithms 
Residual error in pressure as defined in Appendix J 
Stands for general acoustic scalar field variable 
Duct height 

One of the two imaginary square roots of —1 .0 

Unitvector in duct axial direction 

Number of duct segments 

Unitvector in duct transverse direction 

Wavenumber k = io/c 

Total liner length 

Mach number of mean flow 

Number of modes in eigenfunction expansion of acoustic field in jth 
duct segment 

Acoustic pressure 

Time 

PRECEDMG PAGE BLANK NOT FILMED 


ix 


NOMENCLATURE (CQNT'D) 


Axial acoustic velocity component 

Vector acoustic velocity 

Transverse acoustic velocity component 

Transverse duct coordinate 

Axial duct coordinate 

Liner specific impedance 


Real part of modal eigenvalue 

Imaginary part of modal eigenvalue 

Constant defined in Appendix G 

Transverse component of acoustic particle displacement 

Axial velocity component contributed to right moving wave system 
by nth eigenfunction (Appendix F) 

Ratio of duct height to wavelength of sound, ^ = H/X 

Wavelength of sound, ^~2 7T c/U) 

Eigenvalue 

Dimensionless propagation constant, V= XfU. 

Quantity defined in Appendix A 

Axial velocity component contributed to left moving wave system by 
nth eigenfunction (Appendix F} 

Ratio of the length of the circumference of a circle to the diameter 

Axial acoustic energy flux 

Mass density of fliud 

Constant defined in Appendix G 
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Denotes left moving wave 
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THEORETICAL DEVELOPMENT 


Problem Formulation 

An analytical technique has been developed for studying the propagation of sound waves in flow 
ducts with segmented liners. The configuration investigated consists of a uniform geometry 
two-dimensional duct containing uniform mean flow with the duct linings and the acoustic 
excitation symmetric about the duct centerline. The acoustically lined segments were placed in an 
infinite hard-walled duct and the reflection and transmission characteristics resulting from incidence 
of a specified combination of duct modes were investigated. Acoustic fields for ducts with 
nonuniform boundary conditions have been described by eigenfunction expansions {refs. 7-9), 
perturbation techniques (refs. 19-21), finite difference formulations {refs. 22-23), and variational 
and weighted-residual approaches {refs. 10, 11, 24, 25). The geometry of the present problem 
makes it amenable to efficient solution by the mode matching or eigenfunction expansion method. 
The problem was, therefore, formulated as a multiregion coupled boundary value problem. Since 
the boundaries of the duct correspond to coordinate surfaces of a rectangular coordinate system for 
each region, the acoustic fields may be expressed In terms of eigenfunctions which explicitly satisfy 
the Helmholtz equation and the wail boundary conditions. In addition, to ensure continuity of the 
field, solutions in each region must have coinciding boundary data at the junction between regions. 

Acoustic Field in Duct Segment 

Duct Eigenmodes 

The acoustic field in a two-dimensional flow duct is governed by the convected wave equation and 
by the axial and transverse components of the momentum equation: 
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The primed symbols used to denote the physical quantities in these equations are functions of y, z 
and t. The perturbation velocity V', a vectorfield, is defined in terms of its axial and tranverse com- 
ponents. 


V' = v'j + w' k 




(4) 
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where the axial component is expressed in terms of axial particle displacement 
The axial wave propagation and exponential time dependence are characterized by: 


iwt-iXz 


g'(y,z,t) = 

where g' stands for fields p', v', v/, or and unprimed quantities are functions of y alone. 
Use of equation (5) in (1) to (3) yields the following set of reduced governing equations: 


+ k2[(l^ i/M) 2 - p = 0 

0y2 

w = Xp/[pck(1 — V M)] 

V = i4^ / [pckd - pM)] 

dy 

c = _!v /tck(l - i^MH 


The boundary condition at a normally-reacting acoustically lined wall is derived by imposing 
continuity of acoustic particle displacement normal to the wall (ref. 26). This assumption results in 
the expression: 

5P — o 

= - ikA(1 - at wall (7) 

0y 

where A is the specific acoustic admittance of the wall. For the symmetric duct configuration the 
transverse variation of the pressure field is given by the symmetric eigenmodes which are expressed 
in terms of eigenvalues 


Pp = cos ( p^y) 


The eigenvalues p ^ satisfy: 


Pj^tan( PpYl - ikA t V - i^nWl)^ = O 


at the wall, and are related to the modal propagation constants 


P 2 = k^ 1(1 - - p^^] 
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For the hard-wall duct (A = 0.0) the eigenvalues are 

An = 2n7T/H 


{H being the duct height). 


Table 1 contains a summary of the eigenvalue equations and definition of the propagation constants 
for the hard-walled and soft-walled ducts with and without mean flow. Rationale for 
these definitions is given in A,ppendtx A. 

Solutions to the duct eigenvalue equation (8) have been obtained by three different numerical 
techniques: 

• Newton- Raphson iteration {refs. 27-28) 

• Numerical integration of a differential equation which has eigenvalue equation as a 
solution (refs. 4-5) (see Appendix B) 

• Rayleigh-Ritz technique (Appendix C) 

Multiple eigenvalues associated with the Cremer optimum impedance values were determined using 
the method of Tester (ref. 3) (see Appendix D). Zorumski and Mason (ref. 4) have shown that tiie 
conventional duct eigenfunction expansions are no longer valid at this Cremer optimum. To avoid 
having to alter the form of the eigenfunction expansions used in these algorithms, exact Cremer 
impedance values were never used in computing segmented liner performance. It was found that by 
choosing an impedance several percent off the exact Cremer value, adequate eigenvalue separation 
could be obtained to ensure validity of the conventional expansion with little loss in accuracy of the 
performance predictions. 

Eigenfunction Expansions 

The acoustic fields in each of the duct segments are represented by eigenfunction expansions in 
terms of right moving (incident) and left moving (reflected) modes. The axial variation of these 
modes is defined so that they are piiased with respect to the point of generation; right moving 
modes have zero phase at the left end, and left moving or reflected modes have zero phase at the 
right-hand end of the segment. 

Modal expansions for the incident hard-wall, a typical soft-wall, and the transmitting (semi-infinite) 
hard-walled segments, respectively (Figure 1), are given fay: 


^o“' -TX°z "i^' 

‘^os(2n7ry/H)e * **-S^b°cos(2n7r/H)e ' ^ 


Nj-l 


.] Nj-1 ^] 

. i .-^i _-iX(^(Zj~Lj) 


P;=!Zj 3^cos{/f^y)G * ^ j+2 bjcos(/f^y)6 
J n=o n=o n 


(9) 


( 10 ) 
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TABLE 1 

SUMMARY OF EIGENVALUE EQUATIONS AND PROPAGATION CONSTANTS, SYMMETRIC DUCT MODES 
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SOFT-WALLS 


An ~ 2 n7T /H 
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(11) 


The a° are the known modal coefficients of the incident sound fieid. Expressions for the axial 
acoustic particle velocity in each region are: 


-*“o 

,-iX„ 2 bn A, 
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Matching of Segment Acoustic Fields 

Continuity to first order of mass and momentum between adjacent duct segments can be ensured 
by Imposing continuity of acoustic pressure and axial acoustic particle velocity, respectively. The 
errors in pressure, pi — pi and velocity, wi — wi \ at the duct cross section between segments 
j and I + 1 are minimized by use of the method of weighted residuals. By multiplying the pressure 
residual by the left moving modes in segment] and the velocity residual by the right moving modes 
in segment J + 1, and requiring that the errors be orthogonal to each mode of the expansion, the 
following matching equations are obtained: 



cos(pniY)dy-0, m=0,(1), Nj_“1, 


and 


(151 


-H/2 . 

I (w^-w^'*‘^lcos{/ij^^y)dy= 0 , m= 0 , (l), Nj^.j -1 
-H/2 


(16) 
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This matching procedure involves the tacit assumption that the eigenfunctions form a complete set 
Although completeness can be proved for the no-mean flow modal functions and for the mean flow 
hard-wall duct functions, proof of the completeness of the mean-flow soft-wail eigenfunctions is not 
available. (If they are not complete, the whole concept of eigenfunction expansions becomes 
invalid; thus, they are assumed to be complete,) 

Closed form expressions for the duct matching equations (15) and (16) are given in Appendix E for 
various interfaces in a two-dimensional symmetric duct. These equations are summarized in matrix 
form on Figures 2 and 3 for the case of a semi-infinite hard-wall termination and for an arbitrary 
nonuniform termination, respectively. Development of the arbitrary termination equations is given 
in Appendix F. It was necessary to allow for this case in order to model the termination 
characteristics of the test flow duct facility. The duct termination was designed to minimize 
acoustic reflections for a lower frequency regime, but it was not anechoic at the frequencies used in 
the present investigation. 

Energy Flux Expressions 

Several energy flux expressions are available in the literature for describing acoustic energy 
transmission in ducts with parallel mean flow fields. As noted by Goldstein (ref. 29), the acoustic 
energy flux and acoustic intensity cannot be uniquely defined. There are thus several expressions 
available,. Two expressions were selected for use in this study; a) the form presented by Morfey (ref. 
30), and b) that suggested by Ryshov and Shefter (ref. 31). In terms of the problem variables, these 
are: 



(17) 


where H is the axial acoustic energy flux and the asterisk denotes the complex conjugate. 

A measure of the accuracy of acoustic field matching may be obtained by computing the axial 
energy flux on either side of a junction between duct segments. Both flux expressions have been 
found to provide an accurate measure of field matching, although generally the values computed 
using the Morfey expression match more closely. Unless otherwise noted, all liner performance 
values (sound power attenuations) quoted in this report are based on the Morfey expressions for 
energy flux. 

Appendix G contains energy flux expressions for hard-wall incident, hard-wall transmitting, and 
soft-wall segments. The incident energy flux contains cross product terms due to interaction of 
incident and reflected modes. These terms can contribute to incident energy even if all modes are 
cut off. It is also significant that, in general, the energy dissipation in a soft-walled segment cannot 
be inferred from modal attenuation rates. 
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Figure 3. Matrix Form of Duct Acoustic Field Matching Equations. Nonuniform Termination. 


Semi-Infinite and infinite Duct Liners 


Two additional duct lining models were developed to provide comparisons with published data for 
checkout of the mode matching equations and for comparison of finite and infinite duct theories. 
The fii^ involves the Junction between a semi-infinite hard-walled and a semi-infinite soft-walled 
duct.This case, which allows the estimation of single interface effects and for frequencies below sec- 
ond mode cut-on, providesa comparison with closed form solutions available in the literature (ref. 32). 

The matching equations for the semi-infinite case are: 
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The second model is for an infinite soft-wall (Rice^s model, ref. 33) in which the attenuation 
characteristics of specified initial transverse pressure distributions are investigated. For an infinite 
soft-walled duct with an initial pressure distribution defined by: 



0 

a 

n 


cos(2n7Ty/H) 


(19) 


at z = 0 
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the Fourier decomposition of the soft-wall pressure field is: 


u 

r^\ 


N,-1: 

P =E a 
n:=0 



where 


\)A 


N,-1 

»'=E = 



n=0 




sin{^^H/2) 






In either case, the axial energy flux in the soft-wali duct is: 
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The reference energy flux is given by equation (G1) of Appendix G (with the bp^ 0) for the 
semi-infinite case, but by equation (22) with z= 0 for the infinite duct case. 

Analytical Results 

Computational investigations were performed to identify attenuation mechanisms and define 
optimum segmented liner configurations. Because of the large number of parameters involved, a 
multivariable optimization technique was used to determine segmented configurations that 
produced maximum attenuation for specified modal inputs. The conjugate gradient optimization 
technique was used, with the total sound power attenuation as a cost function, and various 
constrained segment lengths and resistive and reactive impedances as independent variables. 
Excellent results were obtained, but it should be pointed out that the optimum liners referred to 
here are generally local optimums but not necessarily global optimums. Also, the local optimum 
obtained in any given solution may vary depending on the starting values and the optimization 
variable step sizes. 
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In all of these studies the acoustic source was specified by a set of incident modai pressure 
amplitudes in the left-hand hard-walled segment To ensure adequate representation of acoustic 
fields, ten terms were used in the eigenfunction expansions [equations (9} through (14j] . Pressure 
and particle velocity profiles were calculated on both sides of each interface to verify that the total 
acoustic fiejds were continuous. Sound power levels were also calculated for both sides of each 
interface and found to rhatch to within a few percent for snail to moderate values of attenuatlbn. 
Soft-wail eigenvalues, modal expansion coefficients, and sound power levels for single and multiple 
segment liners were checked and found to agree with corresponding results obtained by J. F. Unruh, 
who used the mode matching melJiod to study single segment liners at very low frequencies (ref. 
34}. ^ 

Single segment optimum liners were used as baseline configurations against which to measure the 
performance of the multisegment configurations. For the lengths and frequencies under 
consideration, the optimum single-segment liners were identical to those which would be obtained 
using the infinite lined duct theory of Rice (ref. 5), as pointed out in Appendix H. The finite and 
infinite theories produce conflicting results for optirnum single-segment liners only for cases in 
which the excitation frequency is quite close to a cut-on frequency for a low order mode or when 
the segment length is very small compared with the wavelength of the incident sound. For the latter 
cases more sound can be reflected by low resistance liners using length resonance effects than can be 
attenuated by high resistance disipative liners. 

Two-Segment Optimum Liners 

Figure 4 shows the results of a two-segment plane wave constrained optimization run. For this case 
the total lining length and the second segment impedance were held constant while the first 
segment length and impedance were allowed to vary. The resistance of the first segment was 
constrained so as to maintain a reactive front liner. Total attenuation for the single-segment 
optimum for this case was 17.9 dB;the attenuation provided by the two-segment liner was 28.9 dB, 
obviously a substantial improvement. One of the most significant discoveries made during the 
multisegment studies is well illustrated by this case. On Figure 4 the attenuation contributions due 
to reflection from the first interface, absorption in the first segment, and absorption in the second 
segment, respectively, have been plotted separately versus iterations in the optimization procedure. 
As shown, there is effectively no reflection taking place. Increased attenuation for the final 
configuration is due solely to the fact that the first liner has varied the modal content (or, 
equivalently, the transverse pressure distribution] of the sound incident on the second segment. This 
made the second segment more effective. Figures 5 through 10 show the results of an investigation 
into the nature of the modal redistribution caused by the front liner of this two-segment 
configuration. 


Detailed pressure and velocity data at the inlet side of the dissipative liner were computed for both 
the initial /^a" and final "b" two-segment configurations defined on Figure 5, Here "a" refers to the 
configuration chosen as the Initial guess to the optimum liner and "b" to the configuration 
identified by the optimization algorithm as the best liner for attenuation of the incident plane wave. 
In addition, the dissipative liner was placed as a single-segment into a hard-walled duct (Figure 5) 
and subjected to various incident modal combinations to determine the source configuration for 
which the liner produced maximum attenuation. Figures B and 7 illustrate the results obtained, and 
show that optimum attenuation was obtained for the case with the first and second modes incident 
with 48 percent of the incident energy in the first mode, 52 percent in the second, and with the 
second mode phased 100 degrees with respect to the first mode. 
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TWO SEGMENT CONSTRAINED OPTIMIZATION 
PLANE WAVE INPUT, kH ^ 10.0 


INDEPENDENT VARIABLES; L^.Z^ 

CONSTRAINTS: .15H 4 L^ ^1,5H 

REAL(Z.]) :^0.1 

FIXED PARAMETERS: L^ + Lg = 3.0H 

Z2 = 1.581 - 1.3401 
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Figure 4. Determination of Constrained Optimum Two-Segment Plane Wave Liner Using Conjugate 
Gradient Optimization Technique, M = 0.0, 



— -|.Sh[' 2.5H 



INITIAL TWO SEGMENT DUCT, PLANE WAVE INCIDENT, 
= .OS - 1.01, Zj = 1.S81 - 1.341 


b. 


.526H— I ^474H j 


FINAL TWO SEGMENT DUCT, PLANE WAVE INCIDENT, 
Z, = .079 - 1.4051, Zj= 1.581 - 1.341 



I 2.474H 



single segment DUCT, LINER = 2ND SEGMENT OF CONFIGURATION b, 

Z= 1.581 -1.34i, SUBJECT TO INCIDENT MODES FOR WHICH THIS SINGLE 
SEGMENT PRODUCES MAXIMUM ATTENUATION, 47.6% OF INCIDENT 
energy in PLANE wave. 62.4% IN SECOND MODE, SECOND MODE 
PHASED 105° WITH RESPECT TO PLANE WAVE 

Figure 5. Two-Segment Liner Configurations Used in Attenuation Mechanisms Study. 



MODAL PHASE ANGLE, DEGREES 
(2ND MODE WITH RESPECT TO 1ST) 

Figure 6. Performance of Aft Section of Two-Segment Liner as Function of 
Equivalent Hard-Wall Modal Distribution of Incident Sound, M = 0. 
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PERCENT OF INCIDENT ENERGY IN 2ND MODE 

Figure 7, Performance of Aft Section of Two-Segment Liner as Function of 
Equivalent Hard-Wall Modal Distribution of Incident Sound, M - 0. 



Figure 8. Transverse Distribution of Pressure at Entrance to Aft Duct 
Segment for Duct Liner Configurations of Figure 5, M = 0. 
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Figures 8, 9 and 10 show the transverse pressure, axial velocity, and product of the pressure and the 
conjugate of the axial velocity, respectively, at the incident plane of the dissipative liner. These 
quantities, normalized witii respect to their centerline values, are plotted in each figure for cases a, b 
and c. In going from configuration a to b, the optimization algorithm forced the acoustic field 
incident on the liner into a close approximation of the modal content for which the dissipative liner 
performs best. 

Thus, the improved performance is seen to be due to redistribution of the modal content of the 
sound incident on the dissipative liner and not the result of reflection of acoustic energy. (It should 
be pointed out, however, that the method of operation of the low resistance front liner appears to 
be through interaction of right and left moving wave systems; thus, modal reflection in the reactive 
liner is a key to the phenomenon but energy reflection from the reactive liner is not.) The modal 
redistribution mechanism has also been identified by Baumeister (ref. 35). 

Table 2 contains optimum two-segment plane wave liners found by taking various initial guesses and 
constraining only the minimum resistance of the front liner and total length, in each case, the final 
local optimum configuration obtained had a low front liner resistance, even though the initial 
iteration was begun with substantial resistance values. Similar trends have been found by Baumeister 
(ref. 35) and by Quinn (ref, 36) in their investigations. This is gratifying in view of the significant 
differences between their finite difference models and this modal expansion model. (The finite 
difference model involves a finite segmented liner bounded on the incident side by a semi-infinite 
continuation of the initial liner, and on the transmitting side by a pc termination, while the present 
eigenfunction expansion model uses hard -wall terminating ducts on both sides of the lined 
segments.) However, it is important to note that, because of the model difference, these optimum 
configurations will have different impedances and different attenuations for the liner configurations 
in which reactive elements are found. (These configurations include not only the two-segment liner, 
but also the single-segment liners for those cases for which reactive effects are important.) The most 
significant effect of the difference in models will be that, for example, an initial plane wave input 
may be distorted by higher order mode reflection in the model of this report, while the distribution 
will remain plane in the finite difference models since they do not account for the impedance 
discontinuity. 

figure 1 1 depicts the variation with frequency of optimum segment impedances and lengths for a 
two-segment plane wave liner, with L/H = 3.0 and zero i\/>ach number. Specific acoustic resistance 
values were constrained to be greater than 0.048, but the variables were otherwise not restricted. 
The frequency range 6.0 < kH <10.4 is covered. For frequencies greater than kH = 9.6, the 
optimum front liner goes to minimum resistance. Maximum attainable sound power level 
attenuations decrease with increasing frequency. (Although not shown, they are well over 1 00 dB 
for frequencies below 3.0; they drop to 30 to 50 dB in the frequency region investigated.) These are 
definite improvements to be gained in the frequency range studied; by superimposing the single 
segment attenuations on the figure, it is seen that 10 to 15 dB may be achieved by using a 
two-segment configuration. 

Also shown on Figure 11 is a curve of Cremer optimum impedance for the (0, 1) mode pair plotted 
versus reduced frequency. For a reduced frequency of kH near 10, the impedance of the dissipative 
liner, Z 2 approaches that of a Cremer liner. For this case, the front liner apparently converted the 
transverse acoustic field incident on the dissipative liner into a shape similar to the least attenuated 
mode. This behavior does not occur for reduced frequencies kH< 10, and an examination of Table 
2 reveals that local two-segment optimums can be found for kH = 10 that do not approach a Cremer 
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Figure 9. Transverse Distribution of Velocity at Entrance to Aft Duct Segment for Duct Liner 
Configurations of Figure 5, M = 0. 



Figure 10, Transverse Distribution of Product of Pressure and Conjugate of Velocity for 
Duct Liner Configurations of Figure B, M = 0. 
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TABLE 2 

OPTIMUM TWO-SEGMENT LINING CONFIGURATIONS FOR PLANE WAVE, L/H = 3.0, 77 = 1.8, M = 0, 


RUN 

IDENTIFICATION 


Liner 1 


Liner 1 
Liner 2 


INITIAL CONFIGURATION 

IMPEDANCE 

LENGTH 

L/H 

ATTENUATION 

,079-1.40i 
1.581 - 1.34i 

.526 

2.474 

28.9 

1.0 -1.40i 

1.581 - 1.34i 

.526 

2.474 

15.5 

2.0 -1.40i 

1.581 - 1.34j 

.626 

2.474 

13.3 

.079 ■ 1.401 
1.581 - .1i 

.526 

2.474 

11.8 

1.0 -1.40i 

1.581 - .1i 

.526 

2.474 

8.4 


FINAL CONFIGURATION 


LENGTH 

L/H 


.048-1.75i .740 

1.424- 1.1 9i 2.260 


.048-1.70 


.048-1.77i 

1.400-1.211 


.048-1.52i 
1.753- .50i 


.048-1.671 
1.500- .930i 


.829 
2 171 
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: Figure 11. Attenuation, Segment Lengths and Impedance for Optimum 2-Segment Plane Wave Liner as 

J Functions of Reduced Frequency. Total Liner Length = Three Duct Heights. 




value for the dissipative liner impedance. Furthermore, variation of front liner reactance was found 
to be contrary to the results of Baumeista* (ref, 36) which show that the reactance of the front liner 
should also approach the Cremer or least attenuated mode value. 

Three-Segment Optimum Liners 

The analytical studies which were conducted to select lining configurations for the experimental 
tests revealed that the attenuation performance of optimum two-segment lining configurations is 
sensitive to the modal content of the incident wave. Three-segment configurations were studied to 
determine whether an additional segment could reduce this sensitivity to modal content and also 
provide more attenuation than a single segment liner. In conducting the optimization studies, 
optimum single segment plane wave liners for the lengths L-j/H = L 3 /H = .5 and L 2 /H = 2.0 were 
used as Initial values. Final configurations obtained by using these single segment optimums as 
initial values are shown in Table 3. Also shown are configurations that were obtained from other 
arbitrary starting values of liner impedance for a plane wave input. The optimization program 
located many local attenuation peaks that have nearly identical attenuations. Note that for the final 
configurations for cases A-E, the front liner has a high resistance. Only one case could be found for 
the three-segment optimization search which tended toward a low resistance front liner 
configuration. This is different from the trend observed In the case of the two-segment liners. 
However, the three-segment liner also appears to have “modal conditioning" as its operative 
mechanism for improved attenuation. Quinn (ref. 22 ) has obtained similar results; I.e., for his 
three-segment liners, there is no general tendency for the front liner to have low resistance. An 
exception is the case for kH/2 7r = = 1.0 which is the cut-on frequency for the second mode. 

Baumeister (ref. 35) has found that for short liners (low L/H values) the front liner for both two- 
and three-segment optimum liners tend toward low resistance; his data, however, is for r; = 1 .0. As 
noted previously, the single segment liner also tends toward low resistance at ?7 = 1 .0. This repeated 
tendency is apparently due to the ease with which energy is transferred from the plane wave to the 
second mode. It thus appears that near = 1,0 there is some reflection taking place in addition to 
modal redistribution.) 

The final configuration of Case A in Table 3 was used as a starting point in a second optimization 
run. This run was made with different step sizes to determine whether the final configuration of 
Case A could be improved upon. As shown in Table 4, a new three-segment plane wave liner was 
found which out performed the Case A liner. Optimum configurations for other modal input values 
are also shown. 

Table 3 shows that when segment lengths were allowed to vary, several local optimums were found. 
The specific optimized liner parameters, liner lengths and impedances, depended on starting values 
and step sizes used in the search. A distinct global optimum configuration could be found only 
when the segment lengths were held constant. While these local optimum configurations have nearly 
identical total attenuation at the design frequency, they have different sensitivities to Incident 
modal content and different frequency bandwidth characteristics. Because of these properties, it 
may be possible to design a multisegment optimum lining to be relative insensitive to both the 
modal structure of the acoustic field and the tolerances of panel materials. 
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TABLE 3 

LOCAL OPTIMUM THREE-SEGMENT LINING CONFIGURATIONS FOR PLANE WAVE, M = 0. 


RUN 

IDENTIFICATION 


INITIAL CONFIGURATION 


FINAL CONFIGURATION 




Liner 1 
Liner 2 
Liner 3 



Liner 1 
Liner 2 
Liner 3 


Liner 1 
Liner 2 
Liner 3 


Liner 1 
Liner 2 
Liner 3 


Liner 1 
Liner 2 
Liner 3 


IMPEDANCE 

LENGTH 

L/H 

.35- .85i 

0.5 

.87- 1.77i 

2.0 

.35- .85i 

0.5 

2.0 - .85i 

0.5 

.87-1.77i 

2.0 

.35- .85i 

0.5 

.35- .85i 

0.5 

2.0 -1.77i 

2.0 

.35- .S5i 

0.5 

.35- .851 

0.5 

.87 - 1.77i 

2.0 

.05- .85i 

0.5 

.35- .li 

0.5 

.87- 1.77i 

2.0 

.35- .85i 

0.5 

.35- .85i 

0.5 

.87- .It 

2.0 

.35- .85i 

0.5 


ATTENUATION IMPEDANCE 


LENGTH 

L/H 


ATTENUATION 



.048 - 1.419i 
1.869- .744i 
.501 - ,852i 








































TABLE 4 

OPTIMUM THREE-SEGMENT LINING CONFIGURATIONS, M = 0 


PLANE WAVE ' 

AMPLITUDE OF 
FIRST MODE 

1.0 + O.Oi 

EQUAL AMPLITUDE 

1.0 + O.Oi 

LARGE SECOND MODE CONTENT : 

1.0 + O.Oi 

EQUAL ENERGY 

1.0 + 0.0 

EQUAL ENERGY 90° OUT 
OF PHASE 

1.0 + 0.01 


0.0 + O.Oi 


1.0 +0.0I 


5.0 + O.Oi 


1.6034 + O.Oi 


0.0 + 1.6034i 


.982- .977i 
1.571 -2.203i 
.613- .503i 


1.407 - 2.265i 
47 
98 


SEGMENT 

LENGTH 

L/H 


.465 

1.951 

.584 


TOTAL 

ATTENUATION 


.463 

.944 

.594 

53.1 

.536 

.969 

.495 

56.2 


.484 

1.943 

.573 
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EXPERIMENTAL PROGRAM 
Description of Test Facilities 

In the design and subsequent testing of the single and multiple segment lining configurations two 
separate facilities of the acoustic laboratory were utilized. The flow resistance properties of facing 
sheet material were determined in the flow resistance test facility. From the results of these tests, 
facing sheet material was selected for fabrication of the six test panel configurations. The panels 
were then tested in the grazing flow duct facility which provides acoustic data for iining 
configurations In controlled acoustic and flow environmental conditions. 

Flow Resistance Test Facility 

The specific flow resistance test bench, shown on Figure 12, was designed to provide a reliable 
method of determining the steady-state flow resistance of acoustic panel facing sheet materials. The 
components of the specific flow resistance test bench include a device for control of airflow, an 
instrument for measurement of airflow rate, a sealed air duct of known area capable of forcing the 
air through the area of material tested, and an instrument for determining the difference in air 
pressure between the upstream and downstream sides of the material. Output from this 
instrumentation is input to the VARIAN 620L computer which determines values of flow resistance 
in cgs Rayls and calculated values of particle velocity for specified flow rates. 

The accuracy of measurements is estimated to be 5 percent of the measured value in cgs Rayls. This 
estimate is based on the summation of inaccuracies of measuring Instrumentation and of effective 
cross-sectional area of tka specimen tested. Repeatability is within 2 pet •'ent of the determined cgs 
Rayl value. 

Grazing Flow Duct Test Facility 

The grazing flow acoustic impedance test device, schematically shown on Figure 13, consists of an 
air plenum, a sound source, test section, and an aerodynamic/acoustic diffuser. With this device, 
acoustic sound pressure level {SPU and phase data are measured as a function of distance along the 
test section lining length for a selected combination of excitation SPL and frequency and airflow 
rate. 

The air plenum consists of a 121.9-cm {48-inch) diameter cylindrical steel shell lined with acoustic 
insulation to minimize airflow noise. The sound source is provided by an electropneumatic sound 
generator (Figure 14). This generator is coupled to the test section by a fiberglass exponential horn 
{Figure 15) that provides a transition from the circular cross-sectional geometry of the driver to the 
square test section. Sinusoidal noise levels were measured in the test .section for frequencies up to 
6Q00 Hz. These measurements indicate that the sound generator is capable of producing 130 dB for 
Mach numbers up to 0.4. The test section has a 5-cm (2-inch) square cross section with a bellmouth 
inputfromthe plenum (Figure 15). The top of the test section (Figure 16) can accommodate either a 
hard-walled plate or an acoustic liner. The opposite wall has two flush mounted microphones, one 
fixed at the start of the test section and the other on a traversing device (see Figure 17) that moves 
along the 40.6-cm (16-Inch) length of the test section. The bottom and side walls of the test section 
are rigid and acoustically hard. Consequently, the test configuration may be modeled as a 
symmetric duct with total height, H = 10.2-cm (4-inch). Axial flush mounted microphone 
measurements thus correspond to analytical predictions of symmetric duct centerline pressures. 
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Figure 12. Steady-State Flow Resistance Test Bench 
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Figure 13. Schematic Diagram of Grazing Flow Duct Test Facility 
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Figure 14. Installation of the Electropneumatic Sound Generator 


Figure 15. Coupling of the Exponential Horn and Test Section 
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Figure 16. Plan View of Test Section 
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An aerodynamic/acoustic diffuser (Figure 18) is attached to the end of the test section. This 
diffuser with a 3.5-degree diverging angle is designed to allow expansion of the flow field without 
separation. Geometrically, the 5- by 5-cm (2- by 2-inch) cross section expands to a 5- x 34.9-cm (2- 
X 13.75-inch) cross section. The sidewalls are lined with a 50 percent open area perforated sheet. 
The backing space which varies linearly from 0 to 15.2 cm (6 inches) is filled with a bulk material 
acoustic absorber to minimize reflections from the diffuser and provide an anechoic termination. 

A picture of the test operation console is shown on Figure 19. From this console, the acoustic and 
flow environment of the test section are controlled and acoustic and boundary layer data are 
acquired and displayed. The spectral definition, power level, and excitation frequency of the 
excitation system output are controlled by the monitor and control network. The Mach number of 
the test section is regulated by controlling the static pressure of the plenum. 

Design of Test Panel Configurations 


Lining Configurations 

In defining the test panel configurations, it was necessary to design multisegment lining 
configurations which would out-perform single segment linings both for the plane wave and for 
other modal combinations. This behavior was important to the demonstration of the phased lining 
concept since the modal content of the acoustic environment in the grazing flow duct cannot be 
controlled. 

The six test panel configurations are identified and their design conditions specified in Table 5. 
Configurations 1-3 and configurations 4-6 are designed for flow velocities of M = 0.0 and M = 0.4, 
respectively. Configuration 1, a uniform liner, is the single-segment optimum for L/H = 3, and a 
plane wave input for a reduced frequency of kH = 10. For this frequency, (approximately 5K Hz in 
the test section of the grazing flow duct), it is shown on Figure 4 that a low resistance liner placed 
in front of a dissipative liner altered the modal content of sound incident on the dissipative liner, 
thereby making the dissipative liner more effective. This two-segment liner, configuration 3, also 
was selected for the test program. As previously shown, this configuration was a constrained 
optimum; relaxation of the constraints allowed more attenuation to be obtained by subsequent 
two-segment designs. This liner was chosen over the others since a very detailed study of its 
performance had been made. Design parameters for this configuration were determined for M = 0 
by solving equation (23) for facing sheet resistance and backing depth. These parameters are also 
given in Table 5. 

Since the modal content of the grazing flow duct acoustic field was not known at the design 
frequency, the performance of this two-segment liner was calculated for various combinations of 
the first and second modes. The variation with frequency of the attenuation performance of this 
lining configuration was calculated for M = 0.0 by the multisegment duct analysis program. 

Figure 20 depicts the attenuation performance of the two-segment liner. Attenuations for four 
modal combinations are plotted versus excitation reduced frequency. kH. The two-segment liner is 
quite sensitive to the modal content of the sound. The attenuation performances of the 
single-segment optimum lining for a length of 30.5 cm (12 inches) (L/H = 3) are also plotted on tnis 
figure for the same modal combinations. The two-segment liner yields more attenuation than the 
uniform liner only for the case when the plane wave is the only mode present. If there is any 
significant second mode component present in the source field, configuration 3 will not generally 
outperform the single segment liner. 
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Figure 19. Test Operstions Console for Grazing Flow Duct Facility 








TABLE 5 

TEST PANEL CONFIGURATIONS 


CONFIG- DESIGN 

URATION CONDITION PANEL 

I.D. kH M I.D. 


PANEL LENGTH 
cm in 


PANEL 

SPECIFIC RESISTANCE* 
IMPEDWMCE Re PANELS 


FACING SHEET 

FLOW RESISTANCE BACKING 

SHEET SAMPLE Re** DEPTH 
I.D. NO. RAYLS cm in 



18(3 

45 

0*41 

48i6 

41 

0.66 

38(5 

55 

0.36 

18(2 

22 

0.97 

6 

2 

0.61 

48(6 

55 

0.48 

i 


6 58(6 23 


1 




•RESISTANCE VALUES BASED ON pc = 39.5 cgs RAYLS AND REDUCED BY 10 PERCENT TO ACCOUNT FOR CORE AND ADHESIVE BLOCKAGE 
•ESTIMATED VALUES BASED ON FLOW RESISTANCE DATA 
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Obviously, there was a risk of not demonstrating the phased lining concept if the modally sensitive 
two-segment liner were to fail to outperform the single-segment liner due to the presence of the 
second mode In the grazing flow duct. Consequently, the search was continued for a multisegment 
liner that would meet the test objective and would be less sensitive to modal content. A 
three-segment configuration was studied to determine whether the additional segment could reduce 
this sensitivity. In conducting the optimization studies, optimum single-segment liners for the 
lengths L^/H = L 3 /H = .5 and L 2 /H = 2.0, and M = 0 were used as Initial values. 

Figure 21 summarizes the attenuation performance of a three-segment configuration as compared 
with the single-segment optimum. A detailed investigation was conducted to determine the 
mechanism by which the three-segment liner was able to outperform the uniform liner for all four 
modal combinations over the entire frequency range under study. Attenuations for each liner 
segment were calculated for each incident modal combination. For a plane wave input, the 
attenuations of liners 1 and 2 were nearly the same as their optimum single segment attenuation 
values. The attenuation of the third liner is nearly ten times the single segment value of 2.2 dB. 
Thus, in this case, the first two liners accounted for only 5 percent and 29 percent of the total 
attenuation, respectively, but they apparently altered the modal content such that the third 
dissipative liner performed more efficiently {66 percent). Similar behavior can be noted for the case 
with both the first two modes incident at equal pressure amplitudes and equal energies, respectively. 
For these cases, liner 3 accounted for at least 60 percent of the total attenuation. For the case with 
the second mode 90 degrees out of phase with the first (equal energy in both modes), the second 
liner accounted for 53 percent of the total attenuation. It appears that the reason this liner 
performs so well for the various modal inputs is that the function each segment performs, either 
conditioning the modal content or attenuating the principal portion of the sound, depends on 
modal content; the segments change roles for different modal inputs. This three-segment liner which 
w^ found to perform better than the optimum single-segment liner for the plane wave and for 
other selected combinations of duct modes was chosen to be the last IVI = 0 test configuration. 

Configurations 4 through 6 were designed by using the optimization pregram in conjunction with 
the uniform mean flow theory. Plane wave optimum muitisegment lining configurations were 
determined by using optimum single-segment liners for M = 0.4 as initial values. Final parameter 

values are presented in Table 8 . Some of these impedance values are closely approximated by 
1 

applying the j"| -j- factor (ref. 33) to the impedance values of configurations 1 through 3. 

The total length of each panel configuration was 30.48 cm (12.0 inches). Since the moving 
microphone traverses 40.0 cm (15.75 inches), the pressure field can be measured beyond the 
boundaries of the acoustic liners. Configurations 1 through 3 were fabricated with a hard-wall 
section, 10.16 cm (4 inches) in length, in front of the linings to provide measured data across the 
initial hard-wal I/soft-wall interface. This arrangement was modified for configuration? 4 tirrough 6 . 
Hard-walled sections, 5.08 cm (2.0 inches) long, were pliicod in ^ ut of and b 4iind the linings to 
provide measured data across both the initial and final hard-waiiAaoft-waM interfaces. Schematic 
diagrams illustrating the location of the lining configurations and giving the identification number 
for each liner are shown on Figures 22 and 23. Liner segment lengths are presented in Table 5. 

Panel Construction 

The test panels are of conventional sandwich type construction. Arrays of normally reacting 
Helmholtz resonators are provided by small cavities sandwiched between a porous face sheet and 
impervious backing. 
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A fiber metal type of facing sheet was chosen because its flow resistance properties are relatively 
independent of frequency and induced particle velocity. The material used, Brunscoustic plate, 
consists of a stainless steel Brunsmet web that is sinter bonded to a .051-cm (.020-inch) thick 
stainless steel perforated sheet witii a 32 percent open area. All test panels were febrlcated from 
existing stock of 25.4-cm (1 0-inch) by 61-cm {24-inch) Brunscoustic plate. The resonator cavities 
were fabricated from phenolic honeycomb with a .953-cm {3/8-inch) cell size. This core material 
was bonded to tiie facing sheet with a nonflowing roll-on adhesive to minimize adhesive fillet 
blockage of the Brunscoustic plate. 

Selection of Panel Facing Sheets 

To obtain the required flow resistance characteristics, the specific resistance of each panel had to be 
converted into a criterion for the selection of facing sheet material. Panel resistance was calculated 
by multiplying the specific resistance by the value of chartcteristic impedance of air (39.5 cgs 
Rayls). This value was adjusted to account for the blockage {estimated to be 1 0 percent) of the core 
and bonding adhesive. This value was used to select the facing sheet material for panel fabrication. 


Figure 24 depicts steady-state flow resistance data for a typical facing sheet material. Steady-state 
flow resistance values in cgs Rayls are plotted against calculated values of particle velocity. This data 
illustrates the linear behavior of the fiber metel. The small slope of the flow resistance versus 
particle velocity curve allows grazing flow effects to be neglected. Liner impedances are assumed to 
be independent of Mach number. The approximate location of the test areas and the approximate 
locations of the areas that were used for panel fabrication are shown. Each of the other six sheets 
was tested to provide an estimate of Its flow resistance properties. Required values of flow 
resistance and the corresponding measured values for the panel liners are given in Table 5. Figures 
25 and 26 are photographs of the fabricated lining configurations. 


Calculation of Panel Backing Depth 

The backing depth of the resonator cores were explicitly determined from the crossover frequency 
impedance model 


pc pc pc F.H 

w 


- icot QkH){d/H)] 


(23) 


in terms of the face sheet resistance, Rg, characteristic Impedance of air Pc and crossover frequency 
(frequency at which the face sheet reactance equals its resistance) of Fg = 27rfg/c. Unpublished 
impedance tube test data for fiber metal facing sheet materials showed that the reactance of the 
facing sheet varies linearly with frequency. This behavior validates the use of the crossover 
frequency model. Based on these data, a value of crossover frequency fg = 15,000 Hz was used in 
this investigation as a representative value for the Brunscoustic face sheet. 

Analysis of Test Data 

The objective of the test program was to experimentally verify the analytical techniques developed 
to predict liner performance. Duct centerline rms pressure magnitude and phase were calculated 
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Figure 24. Typical Bench Test Steady-State Flow Resistance Data 
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Figure 26. Test Panei Configurations for M = 0.4. 
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using the mathematical model with eitJher a hard-walled or soft-walled panel configuration, and 
these oredictions were correlated with test results. Attenuation performance of each test 
configuration was determined by comparing energy flux before and after the lining. Although 
generally good correlation between theory and experiment was obtained, some portions of the test 
data were not accurately represented by the predictions. To assess the relative level of errors in the 
analytical model and test data, a series of parameter variations and curve fit studies were made for 
the measured data. 

Data Acquisition 

Space correlated measurements of acoustic pressure were made in the grazing flow duct facility for 
each of the single and multisegment lining configurations and a hard-walled test section. These data 
consisted of measuring the amplitude and phase of a moving microphone, relative to a fixed 
microphone, along the test section centerline in one of two test procedures: 

• Frequency Scan 

The separation distance between microphones is fixed and frequency is varied 

• Microphone Traverse 

The separation distance iS varied for a constant discrete frequency 

An extensive test program was conducted and all frequency scan and axial traverse data are 
presented in ref. 37. 

A complete set of acoustic and flow environmentel parameters was measured and recorded for 
each test condition. The temperature of the air supply for the flow duct test section varied with 
ambient conditions. To establish the desired value of reduced frequency, kH, excitation frequency 
was adjusted for each test condition to compensate for temperature changes. To estimate the modal 
content of the acoustic field in the test section for the lining configurations, traverse data was 
measured for a hard-walled test section prior to each of the panel test conditions. 

Frequency scans were made for each panel configuration to experimentally determine the 
frequency of maximum attenuation. Although the panels were designed for a reduced frequency of 
kH = 10 tolerances of the facing sheet material caused some detuning of the liners. To ensure that 
the test panels were tested over a range of frequencies that included their peak attenuations, the 
limits of this range were determined by examination of the frequency scan data. These limits, kH = 
9.6 and 10.16, were then established as the two test frequencies for all panel configurations. 

Determination of Modal Content 

The incident modal content of the duct acoustic field for each test condition was determined by 
analyzing the traverse data for a hard-walled test section. It was assumed that this modal content of 
the incident sound field would not be significantly altered by the presence of the liners. The validity 
of this assumption will be discussed later. The interactions of incident and reflected modes was 
interpreted with the aid of an analytical expression for the quantity pp* whieft is a measure of the 
rms pressure. From the expansion for pressure given by equation (9), pp* is: 
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This expression shows that the rms pressure is a combination of terms which are grouped as eitiier 
interactions of incident modes with themselves and reflected modes with themselves, or Interactions 
between incident and reflected modes. The axial variation of pp* is dependent upon differences of 
propagation constants for the terms that contain the Interactions of incident modes with themselves 
or the interaction of reflected modes with themselves. The terms that contain the contributions of 
th“ ir ions between incident and reflected modes are a function of the sum of the modal 
p constants. An effective wavelength for each of these two types pf Interaction terms can 

L i by dividing 27rby either Real (a^ — or Real (Xp+ This Indicates that 

the axial variation for seif-interaction has a long wavelength while the axial variation for 
cross-interaction has a short wavelength. 


To illustrate the use of equation (24) it is applied to the case of two propagating modes in a 
hard-walled infinite length duct with M = 0 where reflections are neglected. 
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Calculated axial distributions of rms pressure magnitude are shown on Figures 27 through 29 for 
various combinations of fundamental and second Incident modes. The axial variation of pressure 
magnitude is quite sensitive to differences In magnitude and phase of the modal content. Recalling 
that = 2 n 7 r/H for a duct with hard-walls, for a frequency parameter of kH = 10.16 and H = 
10.16 cm (4 inches), the propagation constants become 

Xq = 1.000/cm (2.54/inch) 

X, = .786/cm (2.00/inch) 

The effective wavelength for the interaction of the first and second incident modes is 29.5 cm (11.6 
inches) which correlates well with the oscillations in the rms pressure evident in these plots. 
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Similar reasoning was used to study the interaction between the first and second modes which is 
evident In the measured hard-wall traverses. Measured pressure magnitude (rms) and phase (with 
respect to the fixed microphone) for two typical hard-wailed test conditions are plotted versus 
microphone separation on Figures 30 and 31 for frequency parameters of kH = 9.6 and 10,16 
respectively. In addition to the interaction of the first and second incident modes, interaction 
between incident and reflected modes can be seen to be present for kH = 9.6 as indicated by the 
small amplitude, short wavelength variations superimposed on the larger amplitude long wavelength 
variations. The wavelength of the shortwave Is approximately 3.56/cm (1.4/Inch). The effective 
wavelength for the three possible incident and reflected moda! combinations are: 

Modal Interaction Effective Wavelength 

3.15/cm (1. 24/inch) 

a^ b^ 3.99/cm {1.57/inch) 

3Qb^ 3.51 /cm {1.38/inch) 

These wavelengths indicate that the most probable interaction is between the fundamental incident 
mode and the second reflected mode. Modal reflection is from the impedance discontinuity at the 
junction of the test section and the flow diffuser. 

Modal coefficients for the first and second modes were determined by fitting the traverse data with 
centerline pressure calculated by the multisegment analysis. The analytical predictions are shown on 
Figures 30 and 31 by broken lines. Since the multiple segment analysis does not account for 
reflections at the end of the hard-wall portion of the test section, the analytical predictions do not 
have the small amplitude, short wavelength variation seen in the test data, but overall correlation is 
good. 

The set of modal coefficients determined from the six sets of measured hard-wall data is given in 
Table 6. The magnitudes of these coefficients are nearly the same, while the phase was changed to 
fit the traverse data for the individual tests. 

Figure 32 depicts a typical traverse measurement for a hard-walled test section with uniform flow, 
M = 0.4. Analyses similar to those performed for the M = 0 case were also conducted for uniform 
flow. Very little second mode content was found to be present in the sound field. Therefore, all M = 
0.4 test data were evaluated on the assumption that only the plane wave was incident. 

Analyticai/Experimental Correlations 

Measured liner traverse data were correlated with corresponding analytical predictions to validate 
the improved suppression performance of phased liners. Modal coefficients determined by the hard- 
wall traverse data analyses and liner impedance values based on measured flow resistance data were 
input to the multisegment duct analysis to predict the pressure magnitude and phase along the duct 
centerline for configurations 1 through 6, The use of hard-walled modal coefficients is based on the 
assumption that the presence of a liner does not alter the modal structure of the test section. Liner 
impedances (Tables 7 and 8) for the test frequencies at kH values of 9.6 and 10.16 were calculated 
from measured flow resistance data and the values of characteristic impedance of air, pc, which are 
given for each test condition in Table 9. 
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Figure 31. Comparison of Measured and Predicted Centerline Variations of Pressure Magnitude and 
Phase for a Hard*Walled Test Section, kH = 10.16, M = 0. 
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Figure 32. Measured Centerline Variations of Pressure Magnitude and Phase for a Hard-Wall Test Section, 
kH = 10.16, M = 0.4. 
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TABLE 8 

SOUND POWER LEVEL (PWL) ATTENUATIONS FOR ANALYTICAL/EXPERIMENTAL CORRELATIONS OF 

TEST PANEL CONFIGURATIONS 4 THROUGH 6, M = 0.4. 


LINER CONFIGURATION 

kH 

MODEL CONTENT* 
SECOND 140DE 

SPECIFIC 
IMPEDANCE OF 
EACH PANEL 
SEGMENT 

ATTENUATION 
SOUND POWER 
LEVEL 2idB 

UNIFORM 

CONFIGURATION 4 

9.6 

0.0 + 0.0i 

.640 - 1.036i 

15.1 

10.16 

0.0 + 0.0i 

.640- .927i 

16.5 

TWO-SEGMENT 
CONFIGURATION 5 

9.6 

0.0 + O.Oi 

.046- .7211 
.815- .858i 

21.3 

10.16 

0.0 + O.Oi 

.046- .6391 
.815- .750i 

22.5 

THREE-SEGMENT 
CONFIGURATION 6 

9.6 

i 

0.0 + O.Oi 

.595- .5311 
.778-1.1931 
.416- .245i 

20.4 

10.16 

0.0 + O.Oi 

.595- .4381 
.778-1.0721 
.416- .1581 

26.5 


♦FUNDAMENTAL MODE AMPLITUDE = 1.0+ O.Oi 




CONFIGURATION 
l.D. NO. 


TABLE 9 

ENVIRONMENTAL AND PANEL CHARACTERISTICS 
FOR EXPERIMENTAL/ANALYTICAL CORRELATIONS 


CHARACTERISTIC 
IMPEDANCE OF AIR* 
cgs RAYLS 


.D. NO. 


PANEL 

FLOW RESISTANCE** 
cgs RAYLS 


BACKING DEPTH 



♦BASED ON MEASURED TEMPERATURE AND PRESSURE 
♦♦BASED ON NOMINAL VALUES OF TEST DATA AND INCREASED 10 PERCENT TO ACCOUNT FOR CORE AND 
ADHESIVE BLOCKAGE 
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Measured traverse data for configurations 1 and 4 and the corresponding analytical predictions are 
shown on Figures 33 and 34 for a reduced frequency kH = 9.6 and Figures 35 and 36 for kH = 
10.16. The correlation of pressure magnitude is very good for the kH = 10.16 case. A slight 
difference in phase exists between measured and predicted values near the end of the traverse. The 
measured interactions between incident and reflected waves in the bounding hard-wall sections were 
accurately predicted by the analysis. The correlation for tine 9.6 value of kH is generally poor, 
although the incident/reflected standing wave patterns for the hard-wall sections in front of both 
configurations were accurately predicted. Configuration 4 analytical pressure distributions do not 
agree well with measured data but do exhibit trends similiar to the measured data. The magnitude 
and phase of the predicted axial pressure distribution agree with measured data for configuration 1 
for the first 25.4 cm (10 inches) of traverse distance. However, the measured and analytically 
predicted pressure variations are in poor agreement for the last 12.7 cm {5 inches) of lining length. 
An investigation of possible reasons for this poor agreement is presented later in an evaluation of 
the analytical model. 

Figures 37 through 40 present comparisons of analytical predictions and measured traverse data for 
configurations 3 and 5. The large pressure excursions in the 1Z2-cm (4-inch) hard-wall section for 
configuration 3 were accurately predicted for both test frequencies. These excursions are due to 
interaction of the right moving fundamental mode and the left moving fundamental mode that is 
reflected from the low resistance liner. While the amplitude fluctuations are large, little energy was 
reflected from die interface. Analytical predictions for configuration 5 with an incident plane wave 
do not correlate well with measured data, which would indicate that some second mode is present 
in this test data. In a study to be discussed later, better agreement with the test data was obtained 
by including the second mode in the analytical predictions. 

Experimentai/analytical correlations for configurations 2 and 6 are shown on Figures 41 and 42, 
for kH = 9.6 and 43 and 44 for kH = 10.16. These correlations are quite good considering the com- 
plex structure of standing waves in the three dissipative panels. Reflected wave patterns established 
in the hard-wall sections in front of the lining configurations are predicted very well for each 
corresponding set of hard-wall coefficients. These data illustrate modal conditioning in the first two 
liners and attenuation in the last liner, as was predicted during design of the three-segment 
configuration. The poorest agreement exists for configuration 2 for kH = 9.6. For this case, the 
analysis poorly predicts the behavior of both the phase and magnitude of the measured data in the 
last 12.2 cm (4 inches) of the liner. 

Evaluation of Liner Performance 

While it may be overly optimistic to say that excellent correlation is obtained, it is noted that the 
essential features of the measured data are predicted. This good agreement between theory and 
experimental data indicates that the crossover frequency model is a realistic mathematical model for 
the impedance properties of the test panels and that the modal coefficients determined by analyzing 
traverse data for a hard-walled test section usually provide a reasonable estimate of the modal 
content of the acoustic field for the test section with a soft-wrall. In addition, these correlations 
validate the analysis technique and verify that the previously identified physical mechanisms do 
enable phased linings to provide improved attenuation performance. 

It must be emphasized that the measured axial centerline pressure distribution is not a measure of 
suppression performance of the lining. Sound pressure level attenuation measured for an axial 
traverse indicates duct sound power attenuation only when there is but one soft-wall mode 
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Figure 33. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline Pressure 

Magnitude and Phase for Liner Configuration No. 1, kH = 9.6, M - 0. Incident Modal Coefficients were Determined from 
Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 7 . 
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Figure 34. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Uiner Configuration No* 4, kH = 9.6, M = 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wa!I Data. Panel Impedances are Specified in Table 8. 
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Figure 35. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of 

Centerline Pressure Magnitude and Phase for Liner Configuration No. 1, kH = 10.16, M = 0.0, Incident 
Modal Coefficients were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified 
in Table 7. 
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Figure 36. Analytical-Experinrientai Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 4, kH - 10.16, M « 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 8. 
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Figure 37. Anaiytlcal-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 3, kH = 9.6, M = 0,0. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 7. 
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Figure 38. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 5, kH = 9.6, M = 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wail Data. Panel Impedances are Specified in Table 8. 
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Figure 39. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 3, kH = 10.16, M = 0.0. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 7. 
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Figure 40. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 5, kH = 10.16, M = 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 8. 
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Figure 41 . Analytical-Experimental Correlation- Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 2, kH “ 9.6, M == 0,0. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 7, 


PHASE DIFFERENCE - REFERENCE TO TRAVERSE POSITION 





SPL (dB) TRAVERSE MICROPHONE 




REFERENCE POSITION 


HARD-WALLJLINER NO. 1 
‘ " DF-IF 


AXIAL TRAVERSE POSITION, H ^ 10.16 cm {4 in) 


LINER NO. 2: 
DF-2F 


I < 

LINER NO. 3 I HARD-WALlI-^— 

■ * 


u 

X 


m 


a 


Tl 

"n 

m 

X 

m 

H 

n 

m 


X 


"n 

jn 

X 

m 

Z 

a 

m 

H 

O 


X 

cn 

m 


X 

O 

w 


H 


O 

Z 


o ^ 

§9 

B 5 

5 1 


Figure 42. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 6, kH = 9.6, M = 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel Impedances are Specified in Table 8. 
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Figure 43. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 2 , kH = 10.16, M “ 0.0, incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Pane! Impedances are Specified in Table 7. 
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Figure 44. Analytical-Experimental Correlation. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Liner Configuration No. 6, kH = 10,16, M = 0.4. Incident Modal Coefficients 
were Determined from Analysis of Hard-Wall Data. Panel impedances are Specified in Table 8. 
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propagating. More than one mode can propagate for the test frequencies; consequently, the axial 
standing wave pattern provides a means of correlating analysis and experiment. However, the 
attenuation performance of the lining configurations must be determined by comparing axial energy 
flux before and after the lining. Figure 45 depicts the axial variation of sound power attenuation for 
configuration 4 at a kH = 10.16. A comparison of this distribution of PWL with Figure 36 shows 
that the standing waves present in the SPL distribution would cause attenuation estimates based on 
SPL rather than on PWL to be in error. 

Since limitations of the test facility did not allow measurement of transverse distributions of sound 
pressure level, sound power attenuations were computed for the test conditions. Modal coefficients 
which were determined from the hard-wall data were input to the segmented duct program to 
compute the energy flux values for the liner. If the capability to measure transverse profiles had 
existed, -die energy flux could have been estimated by measuring transverse pressure magnitude and 
phase before and after the liner.* The profile data could have been used to determine modal 
coefficients for the hard-waijed sections or to compute the integral of the rms pressure across the 
duct cross section. (Attenuation performance predicted by this quantity has been found to correlate 
well with the attenuations predicted by the energy flux expression.) 

The sound power attenuations computed for each of the test conditions are presented in Tables 7 
and 8. These calculated values are based on modal coefficients and impedance values obtained from 
test data. Comparisons of the sound power level attenuations of configurations 2 and 3 with 
configuration 1 are shown on Figures 46 and 47. Sound power level attenuations calculated for 
both test conditions are plotted versus reduced frequency, kH. Differences between configuration 1 
{open symbols) and configurations 2 or 3 (closed symbols) illustrate the attenuation improvement 
of phased linings. These plots also show that these attenuations fall near to or within the 
attenuation envelope for the modal combinations used in the design studies (Figures 20 and 21), 

For the zero Mach number tests at kH = 10.1 6, the two- and three-segment liners produced 4.6 and 
6.4 dB more sound power attenuation respectively than the single-segment liner. The two- and 
three-segment liners produced 6.0 and 10.0 dB more attenuation than the single-segment optimum 
for M = 0.4 and kH = 10.16. These attenuation improvements are significant in view of the design 
constraints of the liners. 

Evaluation of Theoretical and Experimental Methods 

Although generally quite good correlation was obtained between theory and experiment, there were 
some portions of the test data which were not accurately represented by predictions. Some of the 
factors which may have contributed to this lack of correlation include the following: 

• The analysis was based on a two-dimensional model while the experimental configuration 
involved a square duct 

^Several measurements of transverse SPL distribution for hard-wall conditions were made in the 
final phases of the program tn an attempt to ftjrther substantiate the modal identification 
procedure. Unfortunately, this limited capability was obtained too late in the program to allow 
application to PWL determination. 
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Figure 45. Calculated Axial Variation of Sound Power Attenuation for Liner Configuration No. 4, kH = 10.16, M = 0.4. 
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Figure 46. An.’tlytical Evaluation uf Configuration 2 Attenuation P:;rformance 
at Test Frequencies. 
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Figure 47. Analytical Evaluation of Configuration 3 Attenuation Performance 
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• Incident modal content varied between tests and with the presence or absence of liners. 

• It was not possible to determine exact values of liner Impedance due to nonuniformity in 
facing sheet resistance, bloeuage effects, and uncertainty regarding the value of crossover 
frequency. 

• The flow duct termination was not anechoic in the frequency regime of the tests. 

Impedance effects were studied by arbitrarily changing the resistance by +1 0 percent for the iVl = 0 
configurations. T.,jse results indicate that while agreement is improved by changes in resistance, in 
some instances initial estimates of resistance were adequate. 

Another trend study shown on previous comparisons was the variation of modal content. This was 
accomplished by trial and error in the early stages of the program and indicated that modal content 
was variant not only from test to test but in the presence of the liner. Insight was gained into the 
actual model content in the presence of the liner as well as into test duct termination effects and 
actual liner impedance values by adjusting these parameters so as to provide a better fit to the 
measured data. Several analytical techniques were developed to automatically and accurately curve 
fit the data. Of those considered, the conjugate gradient optimization method was the most 
successful and will be the only one presented. 

The cost function for application of the conjugate gradient technique is the mean square difference 
between the measured and predicted axial variation of duct centerline pressure (see Appendix J). 
Initial values of incident modal coefficients, and soft-wall segment impedances were input to the 
program, which systematically varied these input parameters to minimize the cost function. 
Modeling of the duct terminating diffuser was accomplished by placing a soft-wall segment of 
uniform Impedance between the semi-infinite hard-wall termination and the test section. Wall 
impedance for this duct segment was also used as an input parameter to the gradient search curve fit 
program. 

Results of this optimum analysis for hard-wail data are shown on Figures 48 and 49. The variables 
for this optimization study were second mode content and termination impedance. Since the 
termination wail Impedance did not change from the starting value of 2 = p c, it was felt that this 
impedance was adequate and it was thus fixed during subsequent analyses. The 
experimental/analytical correlation obtained for the hard-wall data are well within measurement 
accuracies and clearly show the termination impedance trends discussed in previous sections. 

The optimization technique was also applied to configurations 1, 3 and 5 to identify possible causes 
for the poor correlation in some portions of the test data. Application of this technique to the data 
for configuration 1 (Figure 33) was accomplished in two steps. First the data were evaluated by 
varying Incident modal content Little improvement was noted (Figure 50) so next the impedance 
was allowed to vary as well. Results are presented on Figure 51. Here excellent correlation has been 
obtained for both amplitude and phase. Vvnile modal content changed only slightly, a relatively 
large increase in impedance was required to yield the improved correlation. These changes in 
impedance are not within the variations of the facing sheet flow resistance data. Whan attempts 
were made to cut samples from the panel for standing wave impedance tube checks, the face sheet 
separated from the panel, indicating poor bonding of facing to core. Unfortunately, this separation 
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Figi^re 48, Curve Fit to Experimental Data. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Hard-Wall Test Section, kH = 9.6, M = 0.4. incident Modal Coefficients 
and/or Panel Impedance Varied to Give Best Correlation, impedance of Soft-Walled Segment which Models 
Test Section Termination = pc. 
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Figure 49. Curve to Experimental Data. Comparison of Measured and Predicted Axial Distribution of Centerline 
Pressure Magnitude and Phase for Hard-Wailed Test Section, kH = 10.16, M - 0.4. Incident Modal 
Coefficients Varied to Give Best Correlation. Impedance of Soft-Walled Segment Which Models Test 
Section Termination = O c. 
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Figure 50: Curve Fit to Experimental Data. Comparison of Measured and Predicted Axia! Distribution of Centerline Pressure 

Magnitude and Phase for Liner Configuration No. 1, kH = 9.6, Wl = 0.0. Incident Modal Coefficients Varied to Give 
Better Correlation. Panel Impedances are Specified in Table 8. 
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Figure 51 : Curve Fit to Experimental Data. Comparison of Measured and Predicted Axial Distribution of Centerline Pressure 

Magnitude and Phase for Liner Configuration No, 1, kH = 9.6, M = 0.0, Incident Modal Coefficients and Panel 
Impedance Varied to Give Best Correlation. Impedance of Soft-Walled Segment which Models Test Section 
Termination =Pc. 
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precluded running the impedance tube tests, but the improper bonding could allow face sheet 
motion, which would likely cause an increase in both resistance and reactance. This might explain 
the larger Impedance value required for correlation. 

Configuration 3 was also studied by this technique for the kH = 9.6 case (Figure 37). To allow for 
the nonuniformity in flow resistance measured for the dissipative segment (Figure 24), the 
two-segment liner was modeled by a three-segment configuration. The resulting correlation is 
presented In Figure 52 and again is quite good. The impedance values and incident modal content 
are very reasonable and are well within feasible tolerances. Also, the sound power attenuation of 
21.0 dB is the same as tiiat calculated for configuration 3 in Table 7. 

Configuration 5 was studied only for a kH of 10.16. Both incident mode content and segment 
impedance were allowed to vary. The resulting modal content and segment impedance (Table 10) 
are within reasonable estimates of experimental variations. Comparison with traverse data (Figure 
53) shows excellent agreement when the second hard-wall mode is present in the analytical 
predictions. The change in modal amplitude from that measured in the hard-wall case can be 
explained by assuming that the presence of the liner in the test section reflects pressure waves 
toward the source. These waves interact witii the source pressure wave patterns causing a change in 
incident mode content. 

In light of excellent results obtained from the optimized modal contents and Impedances and the 
other previously discussed correlations, it is felt that the two-dimensional duct analyses presented 
herein for the square grazing flow duct are appropriate. 


TABLE 10 

MODAL COEFFICIENTS AND LINER IMPEDANCES FROM AN OPTIMIZATION 
CURVE FIT OF CONFIGURATION 5, kH = 10.16 TRANSVERSE DATA, FIGURE 53 


FUNDAMENTAL MODE ALMPLITUDE 
SECOND MODE AMPLITUDE 
IMPEDANCE FOR LINER NO. 1 


1. + 0.i 


-.18 - .Hi 


.209 - .822i 


IMPEDANCE FOR LINER NO:2 


,713 - .561 i 
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Figure 52: 


Curve Fit to Experimental Data. Comparison of Measured and Predicted Axial Distribution of Centerline Pressure 
Magnitude and Phase for Liner Configuration No. 3, kH = 9.6, M = 0.0. Incident Modal Coefficients and Panel 
Impedance Varied to Give Best Correlation. Impedance of Soft -Walled Segment which Models Test Section 
Termination =Pc. 
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Figure 53: Curve Fit to Experimental Data. Comparison of Measured and Predicted Axial Distribution of Centerline Pressure 

Magnitude and Phase for Liner Configuration No. 5, kH = 10.16, M = OA Incident Modal Coefficients and 
Panel Impedance Varied to Give Best Correlation. Impedance of Soft-Walled Segment which Models Test Section 
Termination - Pc. 









CONCLUDING REMARKS 


It has been shown that multisegment linings can be designed to yield substantially greater 
attenuation than is provided by single-segment liners of the same length. Also, although optimal 
lining configurations vary with the modal content of the Incident sound, thus generally requiring 
accurate specification of the sound source, it was found possible to design a three-segment liner for 
a specific set of conditions which appeared relatively insensitive to modal content. 

It has been found that multisegment liners have "modal conditioning" as their primary mechanism 
for Increased attenuation; the incident sound is redistributed by the front liners into modes which 
are more readily absorbed by the remaining lining segments. These studies have not shown any 
significant reflection of sound energy from the muitlsegment liners (except at a frequency of kH/2 
= 7T, which corresponds to the second hard-wall mode cut-on), even though the semi-infinite 
hard-walled duct on the input side of the segmented liner would readily allow such reflections. 
Other investigators [Baumeister (ref. 35), Quinn (ref. 36), Lester and Posey (ref. 38)1 are in 
agreement with the major conclusion regarding the mechanism by which the multisegment liner 
performance is obtained, even though some of the mathematical models used were quite different 
from those used in this study. 

In addition, a method has been demonstrated for determining duct modal consent for casK with 
only a few modes present The technique involves varying incident and reflected modal coefficients 
to fit measured axial pressure phase and magnitude information. The use of axial trave:se data for 
modal identification in a uniform flow avoids the problems associated with the flow and acoustic 
interference of a transverse probe extended into tiie duct 

Of particular importance is the fact that it was possible to design, build, and demonstrate 
multisegment liners In the presence of such practical constraints as uncertainty regarding the details 
of the sound source modal structure and the requirement that liners be built up from materials in 
stock. (Note that these liners are thus constrained optimums and do not necessarily represent the 
overall optimum liner which might be produced if these constraints could be relaxed.) 

The results of this investigation indicate that the present phased liner theory could be of immediate 
benefit in aircraft applications, particularly to problems of exhaust duct noise suppression. 
However, adaptation to the design of inlet liners will require that several remaining questions be 
answered. These include the following: 

• Can practical phased liners be developed for cases in which large number of modes, 
including spinning modes, are present? 

• How would the presence of sheared mean flow, which affects upstream and downstream 
moving modes differently, change the effectiveness and behavior of multisegment liners? 

• Could the optimization procedure be improved by considering only right moving modes 
in longer, highly dissipative segments, and retaining both right and left moving systems in 
the short and low resistance segments? 

• How well can broadband noise be attenuated? That is, can performance be optimized 
over a finite width frequency band, as well as over model content, liner length, etc.? 
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APPENDIX A 


DEFINITION OF PROPAGATION CONSTANTS 
The propagation constant is defined as 

X = -kM/{1 - m 2) ± k [1 - (^/k)2(1 - m2) ] _ m2) {AD 

The appropriate sign on the square root term must be chosen to designate right and left moving 
modes. 

Soft-Wailed Duct ( H Complex) 

For this case, since ft is complex, define 

[1 - ( (1 - M^)1 ’/==«+ i jS, /?>o (A2) 

Then the axial dependence is obtained by use of (AD, (A2), and equation (5): 

exp [ikM/(l - m 2)+ i « k/(1 - m 2) ± j8k/(1 - m2) ] (A3) 

To ensure attenuation of the soft-wall modes, take the minus sign for right moving modes and the 

plus sign for left moving modes. 

Hard-Walled Duct { fx Real) 

1 

Suppose first that { becomes 

X = -kM/(1 -M2) + ik[(A/k)2n - m2) (A4) 

and the axial dependence is 

exp jikM/d - m 2) ± k ({ ^/k)2 (1 - m 2) _i] -m2) \ (A5) 

Thus, for right and left moving waves, we choose the minus arid plus signs, respectively, in order to 

obtain attenuation of the mode as it propagates away from the source. 

For { // /k)2 <1.0, X is real and the direction of propagation cannot be determined by 
attenuation direction. For this case. 


Pim'i’nixr p; 


' V ''T 


FILMlSSl 
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A=-kM(-M± f )/(l - m2) 


IA6) 


where [1 — ( /r/k)2 (1 — m2)] >M 


Now e 


i £J t - i A 2 




and for wavelike behavior 


Phase j _(] _ o/| ^ + Ml 

Thus the plus sign corresponds to right moving and the minus sign to left moving propagation. 

1 

Finally, consider 1.0 <( /i/k)2<;j . For this case 0<|^ < thus 

A = k(-M + ^ )/{1 - m 2) where 0< M (AS) 


This yields two negative phase speeds for M > 0 


A = - k (M + I )/(1 - M^) 


and two positive phase speeds for M < 0 


A = k( 


M 


)/n -M^) 


(A9) 


(A10) 


Eversman (Reference 39) has shown by energy considerations that the smallest absolute value 
of A corresponds to downstream propagation and the largest to upstream transmission. 

The resulting definitions of flow duct modal propagation constants are summarized in Table 2. 
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APPENDIX B 


NUMERICAL INTEGRATION TECHNIQUE FOR EIGENVALUE COMPUTATION 


This method or similar methods have been used by Rice (ref. 5) and Zorumski and Mason (ref. 4). 


The etgeP'/alue equation for symmetric flow duct modes may be written 


D = H sin ( A H/2)/2 - ikHA cos ( ^ H/2)/2 = 0 


(B1) 


The total differential of D is 


dD = ( dD/dU) df^ +{ dD/ d A] dA = 0 (B2) 

The tqtel differential vanishes because D is identically zero. This equation may be solved for 
d ju /dA: 


d^/dA = [ dD/ d A)nd D/dU) = f{fi ^A) 

= ikH/[ (1 + ikHA) tan { fi H/2) -t /t H] 


(B3) 


This equation may be Integrated from A = 0 to the true softwall value of A using, for example, the 
Runge-Kutta integration procedure 

A(A + <JA) =/^{A) + IK^+ 2 (K 2 +K 3 ) + K 4 J /6 
= 5Af( A) 

Kj = 5Af(A«+K^/2,A+ 5 A/2) (84) 

Kg = <JAf(/^ + K2/2, A+ (J A/2) 

K4 = 5Af(/f + Kg, A+ (JA) 


Note that for A = 0, = 2n tt /H, n is an integer. For n = 0, the differential equation is singular. 

Thus, for the plane wave mode, the integration must be begun at fi {6 A), which may be obtained 
by an asymptotic expansion of (B1) for small admittance or by Newton-Raphson iteration on (B1). 
It should also be noted that, at an optimum Cremer impedance, the quantity d D/df^ vanishes 


(see Appendix D). This method would thus fall if the integration path was to passthrough one of 
these singular points. This possibility could be taken into account in the computational algorithm, 
but would involve considerable program complexity. In experience to date with this program, it has 
been found that the chance of passing through a singular point is rather remote and so the program 
has not been modified. 


APPENDIX C 


RAYLEIGH-RITZ TECHNIQUE FOR EIGENVALUE COMPUTATION 


The functional 
R 


1= p^}dy + iA(1-XMjj)^p^{y=H/2) 


(Cl) 


has stationary conditions 


Pyy + [(1- XMo) 2-X2] P = 0 
Py(y=H/2) + iA(1 -XMo)2p(y=H/2) = 0 
Py(y = 0) = 0 


(C2) 


By expanding p in terms of ftinctions which have the capability of satisfying (C2>, and requiring 
that I be stationary with respect to the undetermined expansion coefficients, a linear algebraic 
eigenvalue problem in can be obtained: 


[l-B] [A] - [U] {a^h {O} 


tC3) 


N-1 


For an expansion p = 2 (n tt y/H), for example, the requirement that 51/ 5 a ^ =0 

. m=o " 

for yC = 0, (1), N ■ 1 leads to the following definitions for the matrices in (C3): 

R 


[A] 


rf 

- J {nmTT 

-n 


^/H^) sin (nTTy/H) sin (myry/H) dy 


R 

— J cos (n tty/H) cos (m tt y/H)dy + iA cos (n tt/2) cos (m7r/2) (C4) 


[-B] 


>4 I 1 

= - 1 


cos (nyr y/H) cos |m yry/H) dy 


Equation (C3) can be solved by any standard eigenvalue package. 
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APPENDIX D 


MULTIPLE EIGENVALUES AND CREMER OPTIMUM ADMITTANCES 


Tester (ref. 3) and Zorumski and Mason (ref. 4) have presented extensive studies of the Cremer 
optimum liner and have shown that the optimum admittance for a given mode pair is that which 
causes the modes to coalesce. The conditions for multiplicity of the eigenvalues are: 


D = AH sin { AH/2) - ikHA cos (AH/2) = 0 
a D/a A = (1+ikHA/2) sin ( AH/2) + AHcos (AH/2)/2 = 0 


(D1) 


or 

sin(AH) + AH = 0 (D2) 

Suppose A = « + ii? , then (D2) becomes 
ctH + sin { C(H) cosh { ^H) = 0 
jS H + cos ( a H) sinh ( ^ H) = 0 


Tester (ref. 3) has shown that for sufficiently high eigenvalues 0. Thus, for ^ >> a 

yS/ce =-cot( aH)tanh{ ySH)— (D4) 

But since tanh ( y5H)->-1.0 for /3 large, (D4) becomes - cot ( ctH) -^0 or 

ctH = m7 T+ tt/ 2, m large integer (D5) 

Note that sin (m tt + tt / 2) =-1'^ and thus (D3) becomes 

mTT + 7 t/ 2 + (“D'^cosh ( ySH) = 0 (D6) 

But cosh > 0 for all ( ySH) thus 

m7T+ 7 r/ 2="(“1)*^[cosh ( jSH)[ (D7) 
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which implies that m must be odd, or m - 2,^+ 1 , thus 


oeH = 2j^7T +37T/2 
and 

jSH = cosh*^ (2 /jt+ 3 7t/2) 

B Logg [ 2 / 3 T+ 3 7t/2+ \/ (2j^7T -3 tt/ 2)^ - i] 
« Logg pjln-i- 3 n/2) + izin + 3 tt/2)] 

for large J 

fiH/2 = (/7r+3 7T/4} + i [ Logg (4/7T+ 37T )]/2 
The optimum admittance is thus 

A = 2i {% + [% - ( M H/2)2] } /kH 


(D8) 


(D9) 


(DIO) 


(DID 


Where + is chosen by the condition that Re (A) > 0, 


APPENDIX E 


INTERFACE MATCHING EQUATIONS 

Equations for the first interface: 

Equal pressure 

Ni-1 

AmC-Z aJ[(1/2)';«nH^(-lAin(Ar^H/2)]/[(^r])^ -(mjr) ] 

n=0 

n,-1 a'L, 

-E bJ[M/2):?rlH='{-irsin(?-’H/2)]e - (mr)"] = -A^a » , 

n=o 

m=0. (1), Nq ~1 

Equal Velocity 

5 {-1 ) "^ sin H/2)] /{(l^n M /k)[(?^ )^“ { nsr)^]! 

Ni-1 

Nf1 ^ 1^1 L (E2) 

e " ’[2cos(//nH/2)cos(7;;^^H/2)][A;;;’tan<;r’ H/2) 

tan H/2)]/ { 1 M/k) [{^T^ f - (^m \ 

Nq-1 

2 0X°[(1/2)7^H^(-l)"sin(?n;H/2)]/)l1^SM/k}[(;2‘n;)2-{n5r)2]}, 

m=0,(1),NTl 

Where 

_ j H, m = 0 

Jh/ 2, m>0 

2 cos(^n ^^2 ) cosl^i^ H/2 ) tan {/Tp H/2)-/?[^ tan H/2)], 

[{?;;)* - 

H/2 + sin { /ip] H/(2^j^) , m=n 
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Equations for intermediate interface: 

Equal Pressure 
Nj-1 


X) aj[2cos(^n H /2 )cos{"?jj^W2)][7n ^ani^n H/2)-J^“ tan(^^ H/2)]/[(^n ^ 
n=0 

-X b_^e " ^[2cos(^n H/2)cos(;3;j"W2l][^i tan{;^^ H 
n=o '' 

— • ^i-l2 

_;^"\an(TrJr’H/2)]/[{;/i)^-<^m ) ^=0. m = 0, (1), N 


Equal Velocity 



N.-1 

1 . . -*-j . ... 

X bj^X^6'^"*”i[2cos{/inH/2)cos{/:i:^H/2)][^n tan(^^H/2) 



tan(^^H/2)]/ 
m 


IE4) 


(1-r„'M/kl [(7n>^-(7ml^][=0. 


m=0, (1), Nj-1 


Where 


2cos(/i:j^ ^H/2 )cos(/a^^H/2) N/2)], 


■'j-1. 


-i-1 


.-^i-1 


-i-i 


-j-1 


B 


mn 


<“n ^ “ l>“m * 

[H/2 + sin {^^*H/2)/(2Xi V nn=n 


m ^ n 
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Equations for final interface: 
Equal Pressure 


N -1 , T^i-1, 

J-2 '^n H-2 2 ^j-2 ir^J-2 -^J-2 ^J-2 -t-J-2 

6 L2cos(^^ H/2)cos{^j^ tan(/Up H/2)-/Zpp, tanIjCZpp H/2)J/ 

rrJ-2)2„ 

Ll/*n ' 

[(2’;[,-2)2-{n!T)^>0, m=0, (1), Nj_2-*1 



Equal Velocity 


N -1 ^J"2, 

J-2 J-2H-J-2 iX^ Lj-2 2 -^I 2 

Xn e [(1/2)7;' V("l)'"sin{Tn“^H/2)]/{(1-Xn M/k)D^ 


J“2 

b 

n=o 


J-2 

n 



[(l/2)^:p| H2(-l)'^sin(^P H/2)]/{(1”Xn M/k)[(^n )^“(mjr)^] 


J-1 ^J-1 

fn Xj^ App|/(1“Xpp M/k)""© , m = o, (1), ^j— 


(E61 


Where 


A 


m ” 


H, m=0 
H/2, m > 0 
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APPENDIX F 


TERMINATION EQUATIONS 


Suppose the duct termination to contain an arbitrary nonuniformity between the last soft-walled 
segment and the final semi-infinite hard-vi^alled duct. Then the acoustic field in the last hard-walled 
uniform segment can be written in terms of the "reflection coefficient" matrix for the 
nonuniformity. Begin, for example, by writing the acoustic field in the nonuniform section in terms 
of right moving and left moving wave systems 

3*n^n (Y'Z) +Db*n ^n (Y-^) 

n n (PI) 

^n 4"n (Y.z)+D b; $ (y,z) 
n n 

where the asterisk denotes the nonuniformity. By equating pressure and velocities on the outlet side 
of the nonuniformity equations for the unknown coefficients can be written as: 

(F2) 

WKl + HKf'fFJV’} 

where the are the coefficients of the transmitted hardwall modes. The may be 

eliminated from equation (F2): 

- WCF]-’[D]]{a*nl + [[B]-[C][F]-'[6]]{b* [ = 

(F3) 

or O2.]k}-fO0]{b;} = )o} 




Next, equate pressures and velocities on the inlet side to obtain: 

[Q] + [R] = [S] {a* } + [T] |b-„( 

Mhn'''^f + M{b„J-2f=[Y]{a*}+[zHb;} 

Using (F3): 

[Q] |a„ J-2} + [KHb„ J-2} = [[T]- [S] [d]-’[6]] |b*„} 


(F4) 


{F5) 


PEECSDING 


^ ^-rr FILMED 
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(F6) 


[w] + M {‘■n-’'"! = [H- M [C’M] l^'n } 

-mK\ 

from which the bp can be eliminated: 

[M-&] W'Kljhn'’'^! +[M-&] [9]'’ W]{bn-'‘"}= I®} (F7) 


{b„J-2}=-[,C]-V]|,„J-2) = [R]ja„J-2( 

(F8) 

where [ft] is the reflection coefficient matrix, which can be derived separately from the segmented 
duct analysis and input directly. Typical methodrfor obtaining [ft] may be found in refs. 7-12 and 
24-25. 







APPENDIX G 


ENERGY FLUX EXPRESSIONS 


Incident in Hard-Wailed Duct 


^ E j [(1 + m2 ) REAL 


-0 

XfnM/k 


^m/k 


1-XmM/k/ ' 1-XmM/k 






1 +m 2 )reAL f— 


1-X“^M/k 




^m/k 


+ REAL 


1-X“ M/k / \ 1-X“ M/k 


+ IVI + M 


+ M + M 


X° /k 


1 -XmM/k 


‘x°^/k 


[»mb^*][(l 




l^'I^M/k 


~Xm/k 




*Xm/k 


^^AL 


X„,/k \* 


l-X^M/k 


+ M.+M 


X° /k 


2 2 \ -*^0 / \ ^ 
' l-X^M/k' \ 1- 


i^m M / 2m7T/kH \ / 2m7T/kH 




Xi^k W 'X"^/k Cm M / ZmTT/kH \ / 2m;r/kH 
' A^ 2 \ / \ i^°^M/k 


■^0 

X^/k 

*;o 

■Xf„M/k 


•t* REAL I a 


,0 b°*r ^ 

'm»m 


">-°mA \'i M , m / Xi'K U "4'K 

^°^M/k / 2 2 \ i:?x‘)„M/k / \ 1 ^°^M/k 


Cm M / 2 m 7 r/kH 


^ \ihX° M/k 




+ JL ( \ / ^m^k Cm j« 

2 I A>/k/\rXmM/k/ 2 


??mM/k ' ^ 


\ . Cm M / 2m;r/kH \l 2m7r/kH 


A>/k/\rx%/k 


Where A =r' ° 

^m ]h/ 2 , m >0 


Jo, m = 0 
f" )H/ 2 ,m >0 
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Transmitted «n Hard-Walled Duct 


reproducibility of iff 
ORIGINAL PAGE IS 


3 2PC^^O / L 


+ M 


' Pxj]]«l/k ' ' Px*^'M/k ' 






’ m=0 


^ ' Sl^/k '^V 


J-1 
-XmM/k 






7 


'm 


2 m 7 r/kH \ / 2 m 7 T/kH \ *ll 

Px'i;WkA 1 ^-j;Wk/Ji 


(G2) 


Transmitted in Soft-Walied Duct 


/7 =:^7.REAL E T.k 


Nj~l N|"l , -UXn-Xm^^i 


a 2PC 


n=0 m=0( ^ ^ 


J ]' 


[M+H+““lFjl + MTiri'][7‘t>n(7„'H/21 

-^i _^j, ir2cos(^nR^2)cos(^Pf,H/2j“| j j* ^ r 

-^ri i*l r— ►]■ ,-*-j -^i* ■ir2cos(^nH/2)cos{;« Jh/2)i 

+wr ][/.„ tanL.„ H/2)-^^tan(,,^H/2^[ J 


m 


j* } 

+a.^b„6 


-i(X„-XnilZj + i>,r^Lj 


(G3) 


m n 


P -_{-^i«'|r 2 coslAnR/ 2 }cos((M;P,H/ 2 j-] fH -«-i 

[M+iwM=ir„-kMr;ri J[ h,; '7 ^ Tv J U„»>><a„h«i 


i, • 1 i *-l. -.-j-.-j.-ir*! — t 

-Am»'’lAiH/ 2 lJ -H=„b„e [M-HWM2|r„-l-MrXJlA„'>"lA„H/2) 


*_j. *|. I r 2 cosl^>/ 2 lcosl#.„;H/ 2 J- 

-A„tan(^„H/2)J[ 2 . . 

I PI ' ' m ^ 


■ d U-. ■ ■■■3 


rtri S':-:] r..^ |,;-J 


Equation (G3) {Cont'd.} 


. ^ ^ f . j, -i(Xn“Xm)2j fp*-]* [yi -^j-#-j 

b'jfe So ^ol‘n^ ® [['■m+M+f rl^ri 

■^i; ,"^l „n , M-^|-*-j*r 25 inl 7 |jH/ 2 )sinlir’H/ 2 l'|r_^i. _ 

-""''.n''^H/2)J+T L - -i" " 


2'n^mJL 


IK 


tanl/t/^H/2) 


"nlH/2j-7r^cot{7T^H/2) 


j im 

e 

n m 


•.. M -:tjitinr^ "°s(;rnH/2)cos(;r;;iH/ 2jir-^j 

M rm + 2 ^n^mJL '‘tV - Ct7j'|2 J L'“n^®"*‘"nH/2) 

■''m'""l/'mH«lJ-^T»o L ' _ ^|.; J[/-„cot(^]H/2J-^ 'co,(7^H/2l] 

n ' A^ rn ' 




+ a'b' e 

m n 


M ,Ti\f T‘rn PL"5 H/2|c^.ir„iH/22-| r i 
2 r mJL Ji^„ 


tan!//^H/2) 


-,.,.aol^„H/2lj^-^<2„a„ L - ^7r - -iV 


• j -4- j 

H/2j-;^^cotl;/nH/2) 


+b'b' e 

n m 


. -i(Xn-Xm)(Zj- Lj) Fr I r2coslTn H/2)cos(;r;;i H/2J- 


M+Fm +~ /'"r 


2 •'m^ n JL J 


1 M-H-i-^j»r2sin(^nH/2)5in(^T^H/2)-lf^i. ^j. ^j ^j *1 1 ) 

i/^n' “ '-“m' J; 


i _ An/k 


<rj = 

1-XiM/k 






where 


APPENDIX H 


COMPARISON OF INFINITE AND FINITE DUCT THEORIES 
FOR SINGLE-SEGMENT OPTIMUM LINERS 


Although E. J. Rice of NASA Lewis (personal communication} has carried out extensive 
investigations of uniform liners optimized for specified input pressure distribution and liner lengths, 
his studies have been based on infinite duct theory. That is, spurious mode generation due to the 
step-function transitions In wall impedance from soft walls to hard walls at the segment ends was 
neglected, and only right-moving waves were considered. As part of the present segmented duct 
study, an investigation was made of single-segment optimum liners including wave reflections from 
the wall impedance discontinuities. 

It was discovered during preliminary data comparisons that excellent correlation existed between 
the finite length single-segment optimums and Rice's plane wave infinite soft-wall duct optimums. 

The data of Figure HI show that the infinite and finite duct results are in good agreement. These 
results raised the question of whether complete modal expansions, including reflected waves, are 
even necessary for dissipative liner combinations which do not include low resistance segments. 

To answer this question, a program was developed to compute attenuations in infinite soft-wail 
ducts for initial pressure patterns corresponding to various combinations of hard-wall duct modes. 
Runs were made with this program and with the semi-infinite hard-wail to soft-wall mode matching 
program to check out the single-segment optimums. Figure H2 shows attenuation contours in the 
impedance plane for the single-segment plane wave optimum as determined using the finite duct 
theory at kH/2 7T = ^7 = 1.6 and L/H = 3.0. In addition, points corresponding to values obtained 
using the infinite duct theory have been identified. These points may be seen to agree quite well 
with the finite theory values. Similar data are shown on Figure H3 which depicts the attenuation 
contours for the case with both the first two modes incident, at equal energy, and in phase. Thus, 
the infinite and finite theories are in good agreement for other modal inputs, as well as for plane 
wave incidence. 

Since these data were all run for kH/2 tt = 1.6, a frequency condition midway between two modal 
cut-on values, it was decided that further runs should be made at kH/2?r = 1.0 and kH/27T= 2.0. 
At these frequencies, which correspond to modal cut-on values, the amplitudes of the spuriously 
generated second and third modes at the hard-wall to soft-wall interface should reach their largest 
values. Figure H4 contains tabulations of attenuation values obtained using the finite, infinite, and 
semi-infinite theories at kH/2 tt = ^ = 1 .0 and L/H = 0.5. Quite poor agreement is evident between 
the semi-infinite and infinite results, and even poorer agreement between the finite and infinite 
cases, indicating that the incident pressure distribution was distorted by spurious generation of the 
second mode at the hard-wall to soft-wall junctions. It should be noted that the finite length 
attenuations indicate that the optimum impedance is tending toward a low resistance (reactive) 
liner, a trend also seen near kH/27r= 1.0 In the single-segment optimizations. This appears to be due 
to the fact that r?= 1.0 corresponds to cut-on of the second mode; at this frequency, it is apparently 
more efficient to reflect energy in the second mode than to di.ssipate it in the plane wave. For plane 
wave incidence at kH/2 7r= 2.0 on a liner of length L/H = 3.0 (Figure H5), there was no essential 
difference between infinite and hard-wali to soft-wall theory predictions. Consequently, it appears 
that there is little tendency for the third mode to be spuriously generated by the lining impedance 
discontinuities. 
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In summary, it appears that infinite or semi-infinite theories may be used for single-segment lining 
design provided that: 

• The liner is highly dissipative (no low resistance segment to cause reflections or standing 
wave patterns) 

• The comolex pressure distribution at the entrance to the segment is accurately 
represented 

• The excitation frequency does not lie near the cut-on frequencies for the low order 
modes 

• The duct length is not extremely small compared with the wavelength of the Incident 
sound 


OPTIMUM 

SINGLE 


SEGMENT 

L/H liners 
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ATTENUATION CONTOURS FINITE THEORY 

■ OPTIMUM FINITE 
• ATTENUATIONS ~ INFINITE THEORY 


Figure H2. Comparison of Infinite and Finite Duct Theories for Single Segment 
Sound Power Level Attenuation Contours, M = 0. 
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ATTENUATION CONTOURS ~ FINITE THEORY 

■ OPTIMUM FINITE 
• ATTENUATIONS ~ INFINITE THEORY 


Figure H3. Comparison of Infinite and Finite Duot Theories for Single Segment 
Sound Power Level Attenuation Contours, M = 0. 
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APPENDIX J 
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tiij 


DETERMINATION OF MODEL COEFFICIENTS 
FROM 

GRAZING FLOW DUCT HARDWALL AXIAL PRESSURE TRAVERSE DATA 


For the hardwall duct configuration, the pressure may be expressed in terms of the following 
expansion: 


N-1 

n=0 


-iX z -iA- 

cos{2n7ry/H)e " + X) ‘^°®^2n7ry/H)e 

n=o 


(z-L) 


{JD 


Where z = 0 at the beginning (reference microphone) measuring point, and z = L is the location of 
the end of the uniform cross section portion of the duct (i.e., the beginning of the duct terminating 
horn). The measured data are comprised of the quantities dB (z) and degree (z) on y = 0 (Figure 
J1): 


dB (z) - 10 Io9io(PwALL<^^ '^WALL ^REF '^REf) 
degree (z) = 180 (^^gp “^^VVALL 


(J2) 


where the traveling probe measures 


P (Z) - 
WALL ~ 


•"WALL^^^ 


-\4> 


WALL 


(z) 


and the reference probe measures 


^RFF " 



e -i 0REF 


(jg) 


(J4) 


Note that (J3) can be written 


‘^WALl(Z) = 10^S<2)/20 


^REF 


i [ TTdeg (Z)/180 - "^REp] 

8 


(J5) 
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DUCT TERMINATION 


L 


Figure J1 . Geometry for Modal identification 


To determine the modal coefficients a^^ and bp, we require that 


p(y “ 0,z) - P^^LL onO<z<Lp{gp 


U6) 


where Lpjgp is the total traverse length. We have chosen to impose {J6) by defining an integral 
squared error: 

•-REF 

e 2 = J [p{y = 0,z) -Pyy^LLz)] [p{y=o,z) U7) 

0 


and requiring that the error be minimized with respect to the coefficients: 

aE^/Sa^ = 0, aE^/ abp^ = 0,m = 0, id, N-1 (J8) 

Alternatively, one could impose 

aE^/a aj = 0 . aE^/ a b^ = o, m= o, id, n-i {J9) 
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Equations (J8) and (J9) are equivalent, producing complex conjugate equation sets for a^^*, 
and a^,, bj^, respectively. We choose to solve the set for a^,, b^,: 




-L) + iX 


m 


z 


z. + i7Tdeg(z)/l80 - i^^g^ 

4 


dz 


(JIO) 


m=0,(l),N-1 



An2 + 


iAllz-L) 


‘m* 


dZ + 


0 



iA (z-L) -J iA *(z-U 

m dz 


, f”^dB{z)/20^iT*{z-U+ i7Tdeg(z)/180- i0-.x 

Pj^^P I 10 e 


m=0, (1), N~1 


(J11) 


The integrals on the right-hand sides of (JIO) and Ull) must be evaluated numerically. Since the 
test data are digitized at equal sampling intervals, the numerical integration was performed using an 
equal interval (Simpson) method. 
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